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The influence of the multi-boson effects on pion multiplicities, single-pion spectra 
and two-pion correlation functions is discussed in terms of an analytically solvable 
model. It is argued that spectacular multi-boson effects are likely to be observed 
only in the rare events containing sufficiently high density fluctuations. The possi- 
bilities of unlike particle correlations for a study of the space-time asymmetries in 
particle production, including the sequence of particle emission, are demonstrated. 



1 Introduction 

The correlations of particles at small relative velocities are widely used to study space- 
time characteristics of the production processes. Particularly, for non-interacting iden- 
tical particles, like photons, this technique is called intensity or particle interferometry. 
In this case the correlations appear solely due to the effect of quantum statistics (QS) 
[1, 2]. Similar effect was first used in astronomy to measure the angular radii of stars 
by studying the dependence of the two-photon coincidence rate on the distance be- 
tween the detectors (HBT effect [3]). In particle physics the QS interference was first 
observed as an enhanced production of the pairs of identical pions with small opening 
angles (GGLP effect [1]). Later on, similar to astronomy, Kopylov and Podgoretsky [2] 
suggested to study the interference effect in terms of the correlation function. 1 

The effect of QS is usually considered in the limit of a low phase-space density such 
that the possible multi-particle effects can be neglected. This approximation seems to 
be justified by present experimental data which does not point to any spectacular multi- 
boson effects neither in single-boson spectra nor in two-boson correlations. These effects 
can however clearly manifest themselves in some rare events (e.g., those with large pion 
multiplicities) or in the eventually overpopulated regions of momentum space; see, e.g., 
[4, 5], references therein and Section 3. 

The particle correlations are also influenced by the effect of particle interaction in 
the final state (FSI) [6, 7]. Thus the effect of the Coulomb interaction dominates the 
correlations of charged particles at very small relative momenta (of the order of the 
inverse Bohr radius of the two-particle system) , respectively suppressing or enhancing 
the production of particles with like or unlike charges. Regarding the effect of the strong 



1 Notc that though both the KP and HBT methods are based on the QS interference, they represent 
just orthogonal measurements [2], The former, being the momentum-energy measurement, yields the 
space-time picture of the source, while the latter does the opposite. In particular, the HBT method 
provides the information about the angular size of a star but, of course, - no information about the 
star radius or its lifetime. 
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FSI, it is quite small for pions, while for nucleons it is often a dominant one due to the 
very large magnitude of the s-wave singlet scattering length of about 20 fm. 

Though the FSI effect complicates the correlation analysis, it is an important source 
of information allowing one to measure the space-time characteristics of the production 
process even with the help of non- identical particles [7, 8]. Moreover, the unlike particle 
correlations, contrary to those of identical particles, are sensitive to the relative space- 
time asymmetries in their production, e.g. - to the relative time delays, thus giving an 
important information not accessible in the standard interferometry measurements [9]. 
In Section 4 we briefly formulate the theory of these correlations and demonstrate the 
possibilities of the corresponding correlation technique. 



2 Formalism 



As usual, we will assume sufficiently small phase-space density of the produced multi- 
particle system, such that the correlation of two particles emitted with a small relative 
velocity in nearby space-time points is influenced by the effects of their mutual QS and 
FSI only. 2 We define the ideal two-particle correlation function R(p\,P2) as a ratio of 
the differential two-particle production cross section to the reference one which would be 
observed in the absence of the effects of QS and FSI. In heavy ion or high energy hadronic 
collisions we can neglect kinematic constraints and most of the dynamical correlations 
and construct the reference distribution by mixing the particles from different events. 

Assuming the momentum dependence of the one-particle emission probabilities in- 
essential when varying the particle 4-momenta p\ and P2 by the amount characteristic 
for the correlation due to QS and FSI (smoothness assumption), i.e. assuming that 
the components of the mean space-time distance between particle emitters are much 
larger than those of the space-time extent of the emitters, we get the well-known result 
of Kopylov and Podgorctsky for identical particles, modified by the substitution of 
the plane wave e WlXl+tp2X2 by the nonsymmetrized Bethe-Salpeter amplitudes in the 
continuous spectrum of the two-particle states tpp^ {xi, X2), where X{ = {U,ri} arc 
the 4-coordinates of the emission points of the two particles and S is their total spin 
[7]. At equal emission times in the two-particle c.m.s. (t* = tl — 1\ = 0) this amplitude 
coincides (up to an unimportant phase factor due to the c.m.s. motion) with a stationary 
solution of the scattering problem ip^_ (r* ) , where k* = p{ = P2 and r* = r* — 
(the minus sign of the vector k* corresponds to the reverse in time direction of the 
emission process). The Bethe-Salpeter amplitude can be usually substituted by this 



2 This assumption may be not valid in the case of low energy heavy ion reactions when the particles 
are produced in a strong Coulomb field of residual nuclei. To deal with this field a quantum adiabatic 
approach can be used [9]. 
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solution (equal time approximation). 3 Then, for nonidentical particles, 

R(puP2) = ^Ps(\^ ) ^*)\ 2 )s- (1) 
s 

Here the averaging is done over the emission points of the two particles in a state with 
total spin S populated with the probability ps, J2s Ps = I- 4 For identical particles, 
the amplitude in Eq. (1) should be properly symmetrized: 

^fltV) [Vf ( k tV) + (-l) S ^, (+) (r*)]/V2. (2) 

Particularly for non-interacting identical particles the characteristic feature of the 
correlation function is the presence of the interference maximum or minimum at small 
|q| changing to a horizontal plateau at sufficiently large |q|. For example, assuming 
that particles are produced independently and that the space and time limitation of the 
production process is effectively described by the one-particle probability of a Gaussian 
form with the corresponding dispersions r$ and Tq , Eqs. (1) and (2) yield: 

R( Pl ,p 2 ) = l + ]T(-l)V;Cxp[-r V 2 - (r 2 +v 2 T Q 2 )q L 2 }. (3) 
s 

Here qr and q^ are the transverse and longitudinal components of the relative 3- 
momentum q with respect to the pair velocity vector v (q T = 2k^ and q L = q a /v = 
27&£, where 7 is the pair Lorentz factor). We see that, due to the relation q = vq = 
vqi , the correlation function at vt > r substantially depends on the direction of the 
vector q even in the case of spherically symmetric spatial form of the production re- 
gion, the interferometric radii squared in the T— and L— directions being rip — Tq and 
r\ = 7q +v 2 Tq . Generally, the directional dependence of the correlation function can be 
used to determine both the characteristic emission time and the form of the production 
region [2]. 

It should be noted that particle correlations at high energies usually measure only a 
small part of the space-time emission volume since, due to substantially limited decay 
momenta of few hundred McV/c, the sources, despite their fast longitudinal motion, 
emit the correlated particles with nearby velocities mainly at nearby points in the 
c.m.s. of particle pair. The dynamical examples are sources-resonances, colour strings 
or hydrodynamical expansion. To get rid of a fast longitudinal motion, sometimes the 
longitudinally comoving system (LCMS) is introduced, in which each pair is emitted 
transverse to the reaction axis. 

3 The equal time approximation is valid on condition [7] \t*\ <C m^lT* 2 for sign(t*) = ±1 respec- 
tively. This condition is usually satisfied for heavy particles like kaons or nucleons. But even for pions, 
the t* = approximation merely leads to a slight overestimation (typically < 5%) of the strong FSI 
effect and, it doesn't influence the leading zero-distance (r* <C |a|) effect of the Coulomb FSI. 

4 For unpolarized particles with spins si and S2 the probability ps = (2S + l)/[(2si + l)(2s2 + l)]. 
Generally, the correlation function is sensitive to particle polarization. For example, if two spin-1/2 
particles are emitted with polarizations Pi and P2 then [7] po = (1— Pi -P2)/4 and pi = (3+Pi -P2)/4. 
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3 Multi-boson effects 

The two-body formalism of Section 2 assumes that the phase-space density / of the 
produced multi-particle system is small compared to unity so that the multi-particle 
correlations can be neglected. The mean phase-space density of the pions of a given type 
at a given momentum p (rapidity y and transverse momentum p t ) can be estimated, 
in the low density limit, as a number of the pions interfering with a pion of momentum 
p and building the Bose- Einstein (BE) enhancement [10]. Using the usual Gaussian 
parametrization for the correlation function in the LCMS: R(pi,p2) = 1 + Aexp(— r 2 x q 2 x — 
r^q 2 , — r 2 q 2 ), and parametrizing the cut-off of transverse momenta by an exponential 
law in transverse mass m t , we can estimate the mean pion phase-space density as 

m x 3/2 cosh V d3n _ ^ V^cxp[-(m t - m)/T] dn 

Wp V <Pp 2 VT(T + m) mi _ dy' [) 

where V = r x r y r z is the interference volume. The mean density is maximal for soft 
pions (pt ~ y ~ 0). At present energies it is typically less than ~ 0.2 thus indicating 
rather small multi-boson effects. These effects can be expected small also in future heavy 
ion exeperiments at RHIC and CERN since presently the LCMS interference volume 
V seems to scale with the density dn/ dy pointing to the freeze-out of the particles at 
a constant phase-space density. Nevertheless, the multi-boson effects can show up in 
certain classes of events or in some regions of momentum space. 

The multi-boson effects can be practically treated provided that we can neglect 
particle interaction in the final state and assume independent particle emission, supple- 
mented by the requirement of a universal single-particle emission function, independent 
of the origin of single-particle sources. The multi-particle emission function then being 
a product of the single-particle ones. The latter allows to calculate the Bose-Einstein 
effect on the original multiplicity distribution or boson spectra with the help of the BE 
weights oj n [4] which can be expressed through the so called cumulants or differential 
cumulants [11]. The calculation of all the cumulants is generally a difficult task. Thus, 
for realistic transport models used to predict particle production in ultra-relativistic 
heavy-ion collisions, the numerical limitations allow to determine only a few lowest or- 
der cumulants (up to about the fourth order) [11]. Fortunately, they are sufficient for 
typical rather moderate pion freeze-out phase-space densities [12]. 

The analytical calculation of the cumulants of all orders is possible in a simple model 
in which the pions are instantaneously emitted according to the Gaussian emission 
function (in Wigner-like phase-space) characterized by the dispersions A 2 and 7q in the 
momentum and ordinary space respectively [4]. 5 

Assuming a Poissonian distribution of the original boson multiplicities with the mean 
multiplicity n, the resulting one, at large n, tends to the BE distribution with the mean 
multiplicity £/(l — £), where £ = 77/(roA+l/2) 3 can be considered as an intrinsic phase 



5 Note that in this model the two-pion correlation function in the low density limit is given by 
Eq. (3) with S = 0, tq = and r 2 . substituted by the dispersion — 1/(4 A 2 ) of the centers of 
the elementary emitters, where 1/(4 A 2 ) represents the minimal possible dispersion of the Wigner-like 
spatial coordinates due to emitter finite sizes. 
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Fig. 1: The intercept of the two-pion correlation functions as a function of the multiplicity n 
and the density parameter £„ lP for several values p — 0, 0.1, 0.2 and 0.4 GeV/c of the mean 
momentum of the two pions. The arrows on the interpolating curves indicate the intercept 
values corresponding to £„ = 3£ = 1.5 (n « 3(n)). 



space parameter. Its maximal value of 1 corresponds to the explosion of the multiplicity 
distribution. Comparing the model phase-space density (/) p =o ~ n/(^/2roA) 3 with the 
experimental estimate of about 0.2, we get for the density parameter £ ss 0.4 — 0.5. It 
was shown [5] that at such densities an approximate ^-scaling of the ratio (n) jr\ takes 
place; (n) becomes close to the asymptotic scaling value of £/(l — £) only at £ > 0.9. 

Regarding the influence of the BE effect on the single-boson spectrum, at sufficiently 
large momenta, when the local density (/) p remains small even at large multiplicities, 
it is dominated by the contribution of the original spectrum, otherwise, at large local 
densities, it is determined by the asymptotic large-density spectrum with the original 
dispersion A 2 substituted by A/(2ro) [5]. Experimentally the effect of BE "condensate" 
was searched for at SPS CERN as a low-pt enhancement, however, with rather uncertain 
results. 

As for the two-boson correlation function, it is well known that for a given multiplic- 
ity n, the intercept R n (0) decreases and the correlation function width increases with 
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the increasing n or decreasing momentum p, both corresponding to the increasing local 
density. In Fig. 1 we show [12] the intercept as a function of the multiplicity n and the 
density parameter £„ iP = exp(— p 2 /2A 2 ), £„ = n/(roA + 1/2) 3 , for several values 
of the mean momentum: p = 0, 0.1, 0.2 and 0.4 GeV/c. As expected, the intercept is 
practically constant at low local densities (small n or high p). As condensate develops, 
the intercept sharply falls down. The sharpness of this drop is however less pronounced 
at higher momenta even if plotted as a function of the local density. Clearly, this lack 
of density scaling is related to a strong decrease of the condensate contribution with the 
increasing momentum. To demonstrate the possibility of the observation of the con- 
densate effect for soft pions produced in certain high multiplicity events not following 
the ordinary proportionality rule between the freeze-out phase space density and pion 
multiplicity in Fig. 1 we indicate by the arrows the intercept values corresponding to 
£„ = 3£=1.5 (n«3<n». 

In the inclusive case corresponding to the original Poissonian multiplicity distribu- 
tion, the correlation function intercept is equal to 2 for any pion momenta or the local 
densities [5] and its width logarithmically increases with the increasing condensate con- 
tribution [12]. 

The results of the considered model should not be taken, however, too literally since: 

a) in contradiction with the experimental indications on a constant freeze-out phase 
space density, in the model there is no correlation between the emission volume and 
pion multiplicity; 

b) the static character of the model can be justified (neglecting the transverse expan- 
sion) in a limited rapidity region only. Thus the pions with a rapidity difference greater 
than about unity have to be considered as originating from different static sources; 

c) since the inclusive spectra at sufficiently high mean multiplicity are dominated 
by the condensate contributions in both extreme cases of very narrow or very wide 
multiplicity distributions of the originally uncorrelated bosons, the intercept of the 
inclusive correlation function is likely to be less than the value of 2 corresponding to 
the Poissonian case; 

f) for identical charged pions, the BE effects are suppressed due to the Coulomb 
repulsion. Since this repulsion is important only in a weakly populated region of very 
small relative momenta, the suppression of the global BE weights w„ is rather small. 
For example, for dj 2 this suppression, being about (ar A 2 ) -1 , is usually less than one 
per mill. The Coulomb distortion of the global multi-boson effects is therefore negligible 
in the rare gas limit. Nevertheless, since the Coulomb repulsion destroys the formation 
of the condensates made up from positive and negative pions in the disjoint phase space 
regions, it can lead to noticeable differences between charged and neutral pions in dense 
systems. Particularly, we can expect a decrease of the chargc-to-neutral multiplicity 
ratio with the increasing phase space density. 

4 Measuring the relative space-time asymmetries 

The correlation function of two nonidentical particles, compared with the identical ones, 
contains a principally new piece of information on the relative space-time asymmetries 
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in particle emission [9]. This is clearly seen in the case of neutral particles when the 
two-particle amplitude ip^_ ^ (r* ) takes on the form 



where the scattered wave </^«(r*), in the considered region of small relative momenta, 
is independent of the directions of the vectors k* and r*. Inserting Eq. (5) into the 
formula (1) for the correlation function, we can see that the latter is sensitive to the 
relative space-time asymmetry due to the odd term <~ sink*r*. Particularly, it allows 
for a measurement of the mean relative delays (t) = (t\ — tz) in particle emission. To 
see this, let us make the Lorentz transformation from the rest frame of the source to 
the c.m.s. of the two particles: r* L = 7(t\l — vt), = Tt- Considering, for simplicity, 
the behavior of the vector r* in the limit \vt\ 3> r, we see that this vector is only 
slightly affected by averaging over the spatial distance r <C \vt\ of the emission points 
in the rest frame of the source so that r* w —jvt. Therefore, the vector r* is nearly 
parallel or antiparallcl to the velocity vector v of the pair, depending on the sign of 
the time difference t = At = t\ — t 2 . The sensitivity to this sign is transferred to the 
correlation function through the odd in k*r* k> — 7k*vi term provided the sign of the 
scalar product k*v is fixed. 

For charged particles there arise additional odd terms due to the confluent hyperge- 
ometrical function F(a, 1, z) = 1 + az/l\ 2 + a(a + l)z 2 /2! 2 + . . ., modifying the plane 
wave in Eq. (5): 



where p = k*r* + k*r*, r\ = (fc*a) _1 , a = {u,ziz 2 e 2 )~ 1 is the Bohr radius of the two- 
particle system taking into account the sign of the interaction (z^e are the particle 
electric charges, u. is their reduced mass), S — argr(l + in) is the Coulomb s-wave 
shift and A c (rj) = 2Trn/{exp(2irn) — 1] is the Coulomb penetration factor. 6 Clearly, at 
a given distance r* , the effect of the odd component in the Coulomb wave function is 
of increasing importance with a decreasing Bohr radius of the particle pair, i.e. for 
particles of greater masses or electric charges. At low energies, the sensitivity of the 
correlation to the odd component can be somewhat modified due to Coulomb interaction 
with the residual charge [9]. 

It is clear that in the case of a dominant time asymmetry, v\{t) \ >• |(7\l}|, a straight- 
forward way to determine the mean time delay (t) is to measure the correlation functions 
i?+(k*v > 0) and i?_(k*v < 0). Depending on the sign of (t), their ratio should 
show a peak or a dip in the region of small k* and approach 1 at large values of k* . 
As the sign of the scalar product k*v is practically equal to that of the difference of 
particle velocities v\ — v 2 (this equality is always valid for particles of equal masses), 

6 This factor substantially deviates from unity only at fc* < 2-7r/|a| (e.g., at k* < 22 MeV/c for 
two protons). Note that for the distances r* > \a\ the confluent hypergeometrical function becomes 
important and compensates the deviation of the Coulomb factor from unity except for the classically 
forbidden region of k* < (|a|r*/2) -1 / 2 , narrowing with the increasing r* . Thus the FSI practically 
vanishes if at least one of the two particles comes from a long lived source (rj, r)' , A, K® , ■ ■ ■)■ 



V'! ( k t ) (r*) = e - k * r *+^(r*), 



(5) 



— ik* r* 



F (-*»?, l,ip) + ^.(r*) 



(6) 
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the sensitivity of the correlation functions R+ and i?_ to the sign of the difference of 
particle emission times has a simple classical explanation. Clearly, the interaction be- 
tween the particles in the case of an earlier emission of the faster particle will be weaker 
compared with the case of its later emission (the interaction time being longer in the 
latter case leading to a stronger correlation). This expectation is in accordance with 
Eqs. (1) and (6) at k* -» 0, (r*) < \a\ and (|^f fe , (r*)|) < 1, when (the arrow indicates 
the limit t>|(i)| > |(r L )|): 

fi + /fi-,l + 2^^1-2MizM 
a a 

The sensitivity of the R + /R- correlation method to the mean relative time shifts 
(introduced ad hoc) was studied [5] for various two-particle systems simulated in Pb+Pb 
collisions at SPS energy using the event generator VENUS 5.14 [13]. The scaling of the 
effect with the space-time asymmetry and with the inverse Bohr radius a, indicated by 
Eq. (7), was clearly illustrated for the K + K~ system (a = —110 fm) and for the like- 
and unlike-sign ttK, np and Kp systems (a = ±249, ±226 and ±84 fm respectively). It 
was concluded that for sufficiently relativistic pairs (71; > 0.5) the R + /R- ratio can be 
sensitive to the shifts in the particle emission times of the order of a few fm/c. Motivated 
by this result the R + /R_ method was recently applied to the K + K~ system simulated 
in a two-phase thermodynamical evolution model and the sensitivity was demonstrated 
to the production of the transient strange quark matter state even if it decays on strong 
interaction time scales [14]. 

The method sensitivity was also studied for AGS and SPS energies using the trans- 
port code RQMD v2.3 [15]. Thus at SPS energy the central Pb + Pb collisions have been 
simulated and the n + K + , ir + p and K + p correlations have been studied [16]. To get 
rid of the effect of a fast longitudinal motion, the study was done in the longitudinally 
co-moving system (LCMS) in which the pair is emitted transverse to the reaction axis 
so that v = v±, tl = Ax and 

Ax* = 7 ± (Ax - v±Ai), Ay* = Ay, Az* = Az. (8) 

The simulated correlation functions R+, R- and their ratios are plotted in Fig. 2. 
We can see that for tt + p and n + K + systems these ratios are less than unity at small 
values of q = k* , while for K + p system the ratio R+/R- practically coincides with 
unity. These results well agree with the mean values of At, Ax and Ax* presented in 
Table 1 ((Ay) ss (Az) w 0). It can be seen from Eqs. (7) and (8) that the absence of 
the effect in the R+/R- ratio for the K + p system is due to practically the complete 
compensation of the space and time asymmetries leading to Ax* w 0. For ir + p system 
the effect is determined mainly by the x-asymmetry. For ir + K + system both the x- 
and time-asymmetries contribute in the same direction, the latter contribution being 
somewhat larger. The separation of the relative time delays from the spatial asymmetry 
is, in principle, possible (see Eq. (8)) by studying the ratio R + /R- in different intervals 
of the pair velocity. 

At AGS energy the Au + Au collisions have been simulated and the ir + p correlations 
have been studied in the projectile fragmentation region where proton directed flow 
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Table 1 Mean values of the relative space-time coordinates in LCMS (in fm) calculated [16] from 
RQMD (v2.3) for 158 A-GeV Pb+Pb central collisions. 
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Fig. 2: Unlike particle correlation functions R+ and R- and their ratios simulated with RQMD 
for mid-rapidity particle pairs 7r + p, n + K + , and K + p in Pb + Pb collisions at SPS energy. 







-o- 








(a) ;i + -p 












-r- * 






(b) ji + - K + 


I.I.I 


I I I _ 


I I I I _ 












(c) K + -p 





Multi-boson effects and space-time asymmetries 



11 



relative to each other both in the longitudinal and in the transverse (flow) directions in 
the reaction plane. It was demonstrated that a modification of the R+/R- method (± 
corresponding now to the signs of the respective components k*) is sufficiently sensitive 
to reveal these shifts [17]. Recently, the shifts predicted by RQMD were confirmed by 
applying the method to the experimental AGS data [18]. 

At low energies, the particles in heavy ion collisions are emitted with the characteris- 
tic emission times of tens to hundreds fm/c so that the observable time shifts should be 
of the same order [9] . In fact the R + / i?_ method has been successfully applied to study 
proton-deuteron correlations in several heavy ion experiments at GANIL [19, 20]. It 
was observed, in agreement with the coalescence model, that deuterons are on average 
emitted earlier than protons. 

5 Conclusion 

Using the analytically solvable Gaussian model, we have discussed the influence of the 
multi-boson effects on boson multiplicities, single-boson spectra and two-boson corre- 
lations, including an approximate scaling behavior of some of their characteristics with 
the phase space density. Though these effects are hardly to be observable in typical 
events of heavy ion collisions in present and perhaps also in future heavy ion exper- 
iments, they can show up in certain classes of events, e.g., in those with high pion 
multiplicities. 

We have shown that unlike particle correlations, compared with those of identi- 
cal particles, contain a principally new piece of information on the relative space-time 
asymmetries in particle emission, thus allowing, in particular, a measurement of the 
mean relative delays in particle emission at time scales as small as 10~ 23 s. To de- 
termine these asymmetries, the unlike particle correlation functions R + and i?_ have 
to be studied separately for positive and negative values of the projection of the rel- 
ative momentum vector in pair c.m.s. on the pair velocity vector or, generally, - on 
any direction of interest. We have discussed here the results of recent studies of these 
correlation functions for a number of two-particle systems simulated with various event 
generators. It was shown that the R+/R- ratio is sufficiently sensitive to the relative 
space-time asymmetries arising due to the formation of the quark-gluon plasma and 
strangeness distillation and even to those expected in the usual dynamical scenarios at 
AGS and SPS energies. As to the detection of the unlike particles with close veloci- 
ties (pi/mi ^2/^2), there is no problem with the two-track resolution since these 
particles, having either different momenta or different charge-to mass ratios, have well 
separated trajectories in the detector magnetic field. For the same reason, however, a 
large momentum acceptance of the detector is required. 

This work was supported by GA AV Czech Republic, Grant No. A1010601 and by 
GA Czech Republic, Grant No. 202/98/1283. 
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